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normal  random  variable  X  may  be  generated  In  terms  of  uniform 
[0,l]  random  variables  u^,  In  the  following  simple  way:  86^ 

of  the  time,  put  X  =  2(u^+U2+Uj-1.5) ,  11^  of  the  time,  put  X  =  1.5(u^+U2-l) , 
and  the  remaining  3'^  of  the  time,  use  a  more  complicated  procedure  so 
that  the  resulting  mixture  Is  correct.  This  method  takes  only  half 
again  as  long  as  the  very  fastest  methods,  is  much  simpler,  and  requires 
very  little  storage  space.  It  is  the  Volkswagen  for  those  who  can’t 
afford  a  Rolls-Royce 

1.  The  Method 

If  you  want  to  generate  a  normal  random  variable  in  a  computer  and 
can  afford  several  hundred  storage  locations,  then  you  should  fashion 
a  procedure  based  on  one  of  the  super  normal  methods,  such  as  the  one 
described  in  [  1  ].  If,  however,  you  are  diffident  about  using  a  super 
program,  then  you  may  wish  to  consider  an  economy  model.  We  will  con¬ 
sider  the  problem  of  producing  a  normal  variable  in  which  ease  of 
programming  and  small  storage  capacity  are  the  primary,  and  running 
time  the  secondary,  considerations.  It  turns  out  that  our  solution  is 
almost  as  fast  as  the  fastest  programs. 

We  want  a  standard  normal  random  variable  X  in  terms  of  inde¬ 
pendent  uniform  [0,l3  random  variables  u,  u^,  U2,...  .  We  seek  a 
simple,  fast  method  for  generating  X  most  of  the  time,  occasionally 
using  a  more  complicated  correcting  procedure.  Our  choice  of  the  easy 

I 


2 


(1) 


method  is  X  =  2(Uj^+U2+Uj-1.5)  •  This  random  variable  is  already  close 
to  normal,  and  we  can  expect  that  ovir  correcting  procedure  will  have  a 
comfortably  small  frequency.  It  turns  out  that  the  greatest  frequency 
with  which  we  can  represent  X  in  this  way  is  .8638,  and  that,  in 
addition,  we  can  carve  out  a  large  portion  of  the  residual  density  by 
putting  X  =  1.5(u^+U2-1)  with  frequency  .1107.  More  specifically, 
let  gj^(x)  and  g^(x)  be  the  densities  of  2(u^-t-U2+Uj-1.5)  and 
l.SCu^+u^-l) ,  respectively: 
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Then  if  f(x)  =  e"^  A/Un) ,  we  have 

f(x)  =  .8638g^(x)+.1107g2(x)+.0228002039gj(x)  +  .002699796lg^(x), 
where  g^  is  the  density  of  the  normal  tail  and  g^  is  the  residual 


density! 
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The  mixture  (1)  provides  us  with  our  method.  We  use  poleu* 
coordinates  to  provide  X  with  density  See,  for  example, 

r  l],  r  2]*  Density  g_  is  drawn  in  figure  It 


GRAPH  OF  gj(x) 


We  may  generate  X  with  density  g^  by  the  rejection  technique  — 
choose  (x,y)  uniformly  from  the  rectangle  and  put  X  =  x  if 
y<g^(x)  or  choose  a  new  (x,y)  if  y>g^(x).  The  efficiency  of 
the  procedure  is  47^,  that  is,  the  probability  that  a  point  (x,y) 
chosen  uniformly  from  the  rectangle  will  lie  under  the  curve  y  = 
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Thus  we  have  the  following  outline  of  a  procedure  for  generating 
a  standard  normal  random  variable  X  in  terms  of  independent  uniform 
[0,1]  random  variables  u^,  u^,  u^,...  : 

1.  With  probability  .8638  put  X  =  2(u^+U2+Uj-1.5) 

(A 

2.  With  probability  .1107  put  X  =  1.5(u^+U2-l) 

5.  With  probability  .0228002039.  fl>rm  pairs  (x,y)  according  to: 

X  =  6u^-3 
y  =  .358U2 

until  y<gj(x) ,  then  put  X  =  x. 

4.  With  probability  .0026997961.  fi)rm  pairs  (x,y)  according  to: 

X  =  v^C[9-2  ln(v^+V2)]/(v^+V2)}^ 

y  =  V2{[9-2  ln(vJ+V2)]/(v^+V2)}^ 

tmtll  either  x  or  y  >3,  then  let  that  one  ^  X.  v^  and 

2  2 

are  uniform  [-1,1],  conditioned  by  ^^+''^2  — 

2,  Comparisons 

A  progreun  based  on  the  above  outline  is  easy  to  write  —  steps  1 
and  2  are  the  important  ones  for  speed,  and  can  be  written  in  machine 
langviage,  while  one  of  the  procedures  such  as  FORTRAN  can  be  used  for 
steps  3  Now  for  a  comparison  of  times  and  storage  space.  On 

the  IBM  7090,  the  above,  written  as  a  function  subprogram  with  all  the 
linkages,  resetting  of  index  registers,  returning  X  in  normeilized 
floating  form,  etc.  takes  about  80  cycles,  whereas  the  super  program  [l],  with 
its  linkages,  etc.  requires  about  50  cycles.  Thus  we  may  produce  normal 
random  variables  in  the  7090  at  the  rate  of  about  6OOO  per  second  with  the 


above  method,  while  the  super  program  in  [l]  will  provide  them  at  the 
rate  of  about  10,000  per  second*  The  Instructions  and  constants  for 
the  short  program  require  about  JOO  words  of  memory  space  in  the  7090, 
whereas  the  super  program  requires  about  IJOO. 

In  the  IBM  1620  decimal  machine ,  the  super  program  takes  about 
20  milliseconds  for  each  normal  variable,  the  above  method,  about 
35  milliseconds,  both  times  being  for  function  subprograms  with  all 
the  necessary  extras  —  linkages,  returning  in  normalized  floating 
point,  two-speed  windshield  wipers,  etc.  The  instructions  and  con¬ 
stants  for  the  short  program  require  about  1000  core  storage  positions 
in  the  1620,  the  super  program,  about  4000. 
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